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Abstract 

We study the Fourier transform of the local density of states (LDOS) in graphene 
in the presence of a single impurity at high magnetic field. We find that the most pro- 
nounced features occur for energies of the STM tip matching the Landau level energies. 
The Fourier transform of the LDOS shows regions of high intensity centered around the 
center and the corners of the Brillouin zone (BZ). The radial intensity dependence of 
these features is determined by the form of the wavefunctions of the electrons in the 
quantum Hall regime. Moreover, some of these regions break rotational symmetry, and 
their angular dependence is determined by the chirality of the graphene electrons. For 
the zeroth Landau level, the ratio between the features at the corners and center of the 
BZ depends on the nature of the disorder: it goes to zero for potential disorder, and 
is finite for hopping disorder. We believe that a comparison between our analysis and 
experiments will help understand the form of the quasiparticle wavefunction, as well as 
the nature of disorder in graphene. 

The Friedel oscillations in the local density of states (LDOS) in graphene in the presence of 
impurity scattering have recently given rise to a lot of theoretical and experimental interest [1, 
2, 3, 4, 5, 6, 7, 8, 9]. The low-energy Friedel oscillations resulting from intra-nodal scattering 1 
decay atypically as 1/r 2 [1, 2], while the oscillations due to inter- nodal scattering decay as 
1/r [4], as expected for two-dimensional systems [3, 10]. This anomalous decay, and in some 
cases the breaking of the rotational symmetry of the high intensity features in the Fourier 
transform scanning tunneling spectra (FTSTS), are manifestations of the chirality of the 
electrons in graphene. These features have already been observed experimentally for epitaxial 
graphene [6]. It is also interesting to note that for exfoliated graphene, the Friedel oscillations 
corresponding to intra-nodal scattering seem to decay as 1/r [8], and not as 1/r 2 , which may 
mean that some internal symmetry-breaking mechanism is at work in this system. 

It appears therefore that the LDOS in the presence of impurity scattering in graphene can 
retrieve important information about the physics of its quasiparticles, specifically about their 
wavefunction. Here we use this to study the wavefunctions of the graphene quasiparticles in the 



x We remind the reader that intra-nodal scattering denotes the scattering of a quasiparticle by an impurity 
such that a small change in its momentum occurs (the quasiparticle does not hop from one Dirac point to 
another), while inter-nodal scattering denotes scattering in which a large change of momentum occurs, such 
that the quasiparticle hops between two Dirac points. Both processes conserve the energy of the quasiparticle. 
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quantum Hall effect (QHE) regime. This issue is of particular importance for understanding 
why the QHE arises in exfoliated, but not in epitaxial graphene, despite the presence of 
Landau levels (LL's) in both their spectra. 

We focus on the regime of strong magnetic field (of order 50T), when the magnetic length 
is of the order of 20 — 30 lattice constants, and when the FTSTS features are easiest to discern. 
For much larger magnetic fields, the corresponding LL energy is too big, and no longer permits 
approximating the quasiparticle dispersion as linear. Furthermore, at huge magnetic fields (of 
order of thousands of T), when the magnetic length and the lattice constant become of the 
same order of magnitude, the intra- nodal and inter- nodal scattering features overlap 2 , which 
makes the FTSTS features harder to interpret. On the other hand, magnetic fields that are 
much smaller (< IT) will give rise to Fourier-space features that are very sharply peaked and 
harder to analyze due to precision limitations. 

For a given value of the magnetic field we calculate the FTSTS spectra when the tip bias 
matches the energy of the Landau levels, focusing in particular on the zeroth and the second 
LL. For energies of the tip situated between two Landau levels, within the approximations 
we use (Landau levels constant throughout the sample, single-impurity scattering, energy 
conservation), the intensity of the FTSTS spectra is negligible. We observe that for the 
zeroth LL, only high-intensity regions corresponding to intra-nodal scattering and to scattering 
between equivalent nodes are present (at the center of the BZ and reciprocal lattice points 
respectively), while no features can be identified at the corners of the BZ corresponding to 
scattering between nonequivalent nodes. The scattering features are rotationally symmetric 
and decay in a Gaussian manner, consistent with the wavefunction of the quasiparticles in the 
zeroth LL of graphene. 

For higher Landau levels, our calculations reveal both intra-nodal and inter-nodal scatter- 
ing features. The intra-nodal ones, as well as the ones corresponding to scattering between 
equivalent nodes are rotationally symmetric. However, the patterns corresponding to scatter- 
ing between nonequivalent nodes are asymmetric. Their asymmetry is a consequence of the 
chirality of the graphene quasiparticles. The radial dependence of these features stems from 
the electronic wavefunction of the Landau levels 3 , and shows intensity minima and maxima. 
The scale of these fluctuations, as well as the scale associated with the Gaussian decay, are 
proportional to the inverse magnetic length. 

The tight-binding Hamiltonian for monolayer graphene is: 



where the operators a\ correspond to creating electrons on the sublattice A and B respec- 
tively, and 



and a is the spacing between two adjacent carbon atoms, which we are setting to 1. 

2 One should note that the distance between these features is given roughly by the inverse of the lattice 
constant, while their size is given by the inverse of the magnetic length. 

3 The intensity is proportional to an integral of a Gaussian and two Hermite polynomials which ends up 
proportional to a Laguerre polynomial. 
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As well known, the energy vanishes at the Dirac points, which are at (see e.g. [11]) 



at — aX 



= £ 5 + / ia i + 



fJLV 

Here £ = ± is the valley index (there are two such nonequivalent points for each elementary 
cell of the reciprocal space), a\ = (2ir/\/3a, 27r/3a), and a* 2 = (—2n/y/3a,2ir/3a), and (/i,^) 
span the family of equivalent Dirac points in the reciprocal space. Note that 

' 01 = ~3~ + 2 / i7r ' K lv' a 2 = g- + 2Z/7T. 

Thus the Hamiltonian can be expanded around the Dirac points k = + q to find 4 

= i?v) (3) 

where v = 3t/(2a). We work with the linearized Hamiltonian, and take into account the 
quasiparticles associated with all the Dirac points in the system. This can be done by adding 
the new indices £, /i, and v to the wavefunction of the quasiparticles, which characterize the 
position of the corresponding Dirac point K^ u . 

In the absence of magnetic field, the eigenfunctions of the above Hamiltonian have been 
extensively studied (see e.g. [13] and references therein). In the presence of a large magnetic 
field (QHE regime), the eigenfunctions have also been determined in Refs. [14, 15] by noting 
that the Hamiltonian reduces to the Hamiltonian of the harmonic oscillator. The diagonal- 
ization of the Hamiltonian for £ = 1, /i — 0, and v = (around the point Kq ) can be done 
by building the eigenfunction [14, 15]: 

V -/L \4>o (y - kl 2 B )) *■ 1 ^V2lW« !/-M»V 

k k,n,a v ' 

where Ib = \/Ti/eB 2Qnm/ y/B\T] is the magnetic length, and <j) n (y)(n = 0,1,2...) = 
e -y / 2l BH n (y) are the eigenfunctions of the one-dimensional harmonic oscillator (H n (y) are 
the usual Hermite polynomials). Also, f = (x,y), the Ck, n ,aS are the annihilation operators 
for quasiparticles in the n+ l'st LL, with wavenumber k along the x direction and band a, and 
Ck-\ is the annihilation operator for a quasiparticle in the zeroth LL. In the new "c"-operator 
basis the Hamiltonian is diagonal, and the Green's functions are: 

G nAa (u) = {c{ yKa {u)c nAa {u)) = 1 , (5) 

U + 10 — & n) k,a 

We generalize this eigenfunction to take into account all the Dirac points (for the first BZ 
this reduces to the wavefunctions described in [15]) and we obtain: 

*(*o= ^>~y*K (6) 



4 Note that here we use a different convention for the definition of the Fourier transform than in [11] which 
yields an opposite sign for f^(k). 
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with 



» A) ~ 2 f2^L\(i + 0My-ki 2 B )J 



/(l + ^ n {y - Hi) - a(l - O0n+l(y - kl%)\ „„ (?) 
^ 2V2L W 1 + O0n+l(y - ^1) + (1 - O0n(?/ - WD,/ M ^ ' 1 J 



The a ^-operators are annihilation operators that beside the wavenumber k, band a and 
LL index n have also the valley indices £, /i and z/. Their unperturbed correlation functions 
do not depend on the valley indices: 

G n,k,a,^) = {4Z^) C n,k,^)) = ~ + ^ _ ^ > ( 8 ) 

and, in the absence of disorder, correlators of operators connecting two different valleys are 
zero. 

We introduce a delta-function impurity localized on an atom belonging for example to the 
A sublattice with an impurity potential 

V = uV A (r = 0)V A (r = 0), (9) 

where \l/(r) is given by Eq.(6). Using Eqs.(6,8) and the Born approximation, we find that in 
the presence of the impurity potential V, the corrections to the correlation functions for the 
c operators are given by: 

-a(l - 0(1 " e)M-kll)<t>n + i(-kT B )\ 

s(< a AU*) = * {oj + iS)(u I iS _ ^ [(i + - OU-kil)M-k'il) 

-a(l - 0(1 - 04>n + i{-kl 2 B )M-k'l B )] 
-a'(l + 0(1 - e)M-kl 2 B )<t>n'+i(-kT B ) 

-a(l - 0(1 + O)0n+l(-^|)^(-^) 

W(i - 0(1 - O^n+^-^D^+iC-fc'zl)] (10) 

We can now use these formulas to compute the corrections to the LDOS due to impurity 
scattering. Given that the expectation value of the density operator is [11, 12]: 

Pa/b(?) = (^ A/B (f)^ A/B (r)) (11) 

with the A and B components of the \l/(r) being given by Eq.(6), we can write the Fourier 
transform of the LDOS as: 

p(q,cu) = J dfe-^[p A {r) + p B {r)e^ 5AB ] (12) 
4 



(the inclusion of the phase factor in the second term is explained in Refs.[ll, 12]). A long but 
straightforward calculation yields for Sp(q) due to impurity scattering: 



/OO l>00 f' £ 1 

dk / dye-^e^'y-^il - £) 2 (1 - O S T -n 2 x 

0o(y - £4)0 o (y - k'l 2 B )M-kl 2 B )M-k'l 2 B )\ kl _ k a , K e K , (13) 
for the zeroth LL (cu — 0), and 

x{0„(y - kl 2 B )(f) n (y - k'l 2 B )[(l + 0(1 + O + e"^ a (l - 0(1 " 01 
-a0 n (y - *4)0 n+1 (y - k'l 2 B )[(l + 0(1 - O - e-** a (l - 0(1 + 01 
-a0 n+1 (y - kl 2 B )Uy ~ k'l 2 B )[(l - 0(1 + O - e^ a (l +0(1 - 01 
+0 n+1 (y - kl 2 B )<p n+1 (y - k'll) [(1 - 0(1 - O + e-^ a (l + 0(1 + 01} * 
x [(1 + 0(1 + O0n(-^|)0n(-^) - a(l + 0(1 " O0n(-A:/|)0„ + i(-A: , /|) 
-a(l - 0(1 + 0<f>n + i(-kl 2 B )U-k'l 2 B ) 

+ (l-O(l-O0n + l(-M|)0 n+1 (-^|)]| fe/=fe _^ } (14) 



for the n + l'st LL and band a (u = aE n ). 

In this calculation we have neglected the contributions coming from a quasiparticle being 
scattered between different LL's, as we work under the simplifying assumptions that the dom- 
inant scattering mechanism is elastic, and that the energy of each LL is constant throughout 
the sample. We have also focused on energies matching the energies of the Landau levels. For 
intermediate energies the intensity of the spectra is greatly reduced in the limit 5 — * 0, as it 
can be seen from the factor of l/[(o> + iS — aE n )(uj + iS — a'E n i)] in Eq. (10). 

For the zeroth LL {u = 0) we compute Sp(q, uj) analytically: 

S P o(q)= £ Spoiq + K^-K-l), (15) 

where 

fOO pOO 

\2 In _/i.l„_i_„ ;, \2 , 



/oo poo 
d ye -'i<iyy e -(y/lB? /2 e -(y/l B +qJb? /"i I ^ e iq y kl% e -(kl B ) 2 /2 e ^(kl B +q x l b ) 2 /2 
-oo J —oo 

0, e ~lU<llW v )/2 _ (16) 



The corresponding spectrum is plotted in Fig. 1. We note that there are regions of high 
intensity corresponding to quasiparticle scattering between two equivalent nodes (at the center 
of the first BZ and all equivalent points related via translation by a reciprocal lattice vector), 
but no high-intensity regions at the corners of the BZ corresponding to scattering between 
nonequivalent nodes. This can be seen directly from Eq. (13) as only the term proportional 
to (1 — 2 (1 — O 2 appears in Sp(q,u) in the zeroth LL, and this term is nonzero only if 
£ = = I. This is related to the fact that in the zeroth LL the electronic wavefunctions 
have only one non-zero component (A or B), depending on the type of node (£ = ±1) on 
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Figure 1: FTSTS spectrum for a monolayer graphene sample with a single delta- function impurity, 
for an energy E = inside the zeroth LL, and Ib/cl = 7. The BZ is indicated by dashed lines. 



which the electron sits. In order for the electron to scatter between two nodes, it needs to 
be able to change the sublattice index during the scattering process. However, for the type 
of impurities we consider (potential disorder), this is not possible. The shape of the observed 
regions is rotationally symmetric, and the intensity decays with the distance from the center 
in a manner characteristic to the decay of the ground state of the harmonic oscillator (e~ x ) . 
For an energy corresponding to the n + l'st LL, the FT of the LDOS is given by: 

Span(q) = £ K%, - K»>»>) , (17) 

where 

/oo 
dye-^{My)My + &J|)[(i + 0(i + O + e-^+^'-^a - 0(1 - 01 
-oo 

-aUv)<t>n + i{y + qJl)[(l + 0(1 - O - e- 2 ^W-')/3 (1 _ 0(1 + ^ 
-acf> n+1 (y)Uy + QJlW " 0(1 + O - e-^W-'J^i + £')] 
+0 n+1 (y)0 n+1 (y + &Z|)[(1 - 0(1 - 0) + e-^+^'-^a + 0(1 + 0)]} 

/oo 
dfar^Ki + 0(i + OM-kiDM-kil + qJl) 
■oo 

-a(l + 0(1 - e)<f>n(-kll)<f>n + i(-kl 2 B + g x Z|) 
-a(l - 0(1 + Otn+ii-klDM-kll + qJl) 

+ (1 - 0(1 - 04>n + i{-kl 2 B )(t>n + i{-kl 2 B + qJ 2 B )} , (18) 

and where we have considered that q is small with respect to 1/a, hence approximating 
e -iq y a ^ i Denoting I mn (q) = e^ iqyV (f) m (y)(j) n (y + q x l B ), we can rewrite the above formula 



as: 

a/C^fl * tt 1 + 2 (i + 2 l«fll 2 + (i - 2 (i - 2 l W+i(g1l 2 

+(i + o 2 (i - riivi®! 2 + (i - o 2 (i + i'f\in + iM\ 2 

-(i - 2 (i + 2 W(?1W«) - (i + 2 (i - O 2 W(g1/; n+1 (g) 

+(i+o 2 (i+n 2 w + i(^^„(^]}- (19) 

We can compute the I mn integrals analytically: 

Imn{q) = V^2 N M\[(a nm qJ B - iq y l B ) /2} N ~ M L N M M \q 2 l 2 B /2)e'^-^ 2 ^ 4 , (20) 

where is a Laguerre polynomial, a nm = sign(n — to), g = |g| is the length of the q vector, 
and M and N denote the smaller and respectively the larger of to and n. 

At this point there are a few observations we can make. The first is that Sp(q, uj) contains 
terms proportional to Laguerre polynomials (which in general characterize the overlap between 
two LL eigenfunctions[16]). A Laguerre polynomial shows a number of zeroes (or "nodes") 
given by the order of the polynomial (e.g. L 1 has 1 node, L 2 has 2 nodes, and so on). As these 
Laguerrre polynomials appear in the dependence of the LDOS Sp(q, uj) on q, we expect Sp(q, uj) 
to also have a "node" structure. Indeed, the scattering between nonequivalent Dirac points 
(£ 7^ £') gives rise to terms that contain only combinations of I n , n +i terms (all proportional to 
a Laguerre polynomial of order n); these combinations will have n nodes. However, scattering 
between equivalent Dirac points £ = £' gives rise to terms that contain both I n and I n+ i, and 
hence contain combinations of Laguerre polynomials of different orders. In general one cannot 
predict the exact number of nodes of such combinations. 

The second observation is that the high-intensity features show, besides the Laguerre 
polynomial dependence, a Gaussian decay on a scale of a few times the inverse magnetic 
length e~ q 1 b/ a . 

The third observation is that the expression in Eq. (20) contains rotationally asymmetric 
phase factors. The terms of the form e iqxqyl B^ 2 cancel in the final expression of p where only 
products of an integral of the type I mn and of a complex conjugate of such an integral appear. 
The other complex phase factors q x ± iq y only appear if m ^ n, and hence, as it can be seen 
from Eq. (19), only for processes involving scattering between nonequivalent nodes (£ ^ £'). 
This type of scattering corresponds to coupling between the LL wavefunctions <p n and <fi n +i- 
For the non-relativistic QHE this coupling can only arise for the transition of a quasiparticle 
between two distinct (energy separated) Landau levels. However, we see that for graphene, due 
to the spinorial structure of the wavefunction, such (0 n , <f> n +i) coupling occurs naturally inside 
the same LL, and appears to be a manifestation of the chirality of the graphene quasiparticles. 

We have already discussed the particular case of the zeroth Landau level n — — 1. To 
illustrate the observations above we focus on the FTSTS spectra in the second Landau level 
(n — 1). This allows us to observe the node structure of the results in more detail than for 
n = (first LL). The I mn integrals are given by: 

I u (q) = v^F(2 - g 2 /|)e-^-^ 2 ^/ 4 

hM = v^fe-^)(4-g 2 /|)e-^-^) 2 ^/ 4 

hM = -v^fe + ^)(4-g 2 /|)e-^-^ 2 ^/ 4 

J 22 (g) = ^[S + qHKqHl-S^e-^-^y 4 (21) 
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Figure 2: a) FTSTS spectrum for a monolayer graphene sample with a single delta- function impurity, 
for an energy equal to the energy of the second LL i.e. uj = 1.4V, and for Is/a = 7, (B = 900T). The 
BZ is indicated by dashed lines, b) Horizontal cut through the spectrum depicted in a) for k y = 
(in arbitrary units). For clarity, the intensity of the features at the corners of the BZ was multiplied 
by a factor of two. 
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In Figs. 2 and 3 we plot the real part of the FT of the LDOS, first for an unrealistically 
large magnetic field, B = 900T (Fig. 2) as well as for a physical B = 50T (Fig. 3). In the 
first case the corresponding energy for the second LL is of the order of 1.4V which makes the 
linear approximation for the spectrum invalid. However we use it to illustrate qualitatively 
our results, as the magnetic length is Ib — 7a, which gives rise to larger and better resolved 
features. For the more physical case of B = 50T, as well as for smaller magnetic fields, these 
features should be observable for energies smaller or equal to 400meV, for which the linear 
approximation is still reasonable [17]. However, at fields of the order of or smaller 50T, the 
magnetic length is of the order of or larger than 26.5a, which gives rise to sharper features, 
therefore harder to resolve in momentum space, as it can be seen in Fig. 3. 




Figure 3: The real part of the FTSTS spectra for a monolayer graphene sample with a single 
delta-function impurity, for an energy equal to the energy of the second LL i.e. oj ~ 400meV, and 
Ib/cl = 26.5 (B = 50T). The intensity of the features at the corners of the BZ was multiplied by a 
factor of two. 



For both values of the magnetic field there exist regions of high intensity corresponding 
to both intra-nodal and inter-nodal scattering which exhibit maxima and minima superposed 
over a Gaussian decay; the distance between the maxima and minima is proportional to 
the inverse magnetic length. The high-intensity regions corresponding to scattering between 
nonequivalent Dirac points (£ ^ £') have only one zero-intensity node (n = 1) in their radial 
dependence, as expected (see Fig. 4 a) and b)). The intensity also goes to zero at the center 
of these regions (q = 0) because of the q x ± iq y factors in Eq.(21). The high-intensity region 
corresponding to scattering inside the same Dirac point (located at the center of the BZ) 
has n + 1 = 2 nodes, while the high-intensity features corresponding to scattering between 
equivalent Dirac points (£ = £') show a structure of minima and maxima similar to that of 
the central feature, but no zero-intensity node. 

Also, as noted above, the features corresponding to scattering between equivalent Dirac 
points are rotationally symmetric (depend only on the magnitude q). However, for scattering 
between nonequivalent Dirac points, p(q) is proportional to (q x ±iq y ) 2 , which breaks rotational 
symmetry. This also happens in the imaginary part of the FTSTS spectrum, depicted in Fig. 4. 
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Figure 4: The imaginary part of the FTSTS spectrum for a monolayer graphene sample with a 
single delta- function impurity, for uj = leV, and Ib/cl = 7 (B = 900T). As before, the intensity of 
the features at the corners of the BZ was multiplied by a factor of two. 



To conclude, we have computed the effect of single- impurity scattering on the Fourier 
transform of the LDOS in the presence of a strong magnetic field. We have found that the 
FTSTS spectra contain high-intensity regions corresponding both to scattering processes in 
which quasiparticles remain at the same Dirac point (intra-nodal), and to scattering processes 
in which quasiparticles hop between different Dirac points (inter-nodal). Both types of pro- 
cesses give rise to features that contain information about the wavefunction of the electrons 
in the quantum Hall state (decay length, maxima, zeroes, minima). While scattering be- 
tween equivalent nodes gives rise to rotationally symmetric features, the features coming from 
scattering between nonequivalent nodes break this symmetry, manifesting the chirality of the 
quasiparticles. A special situation arises for energies inside the zeroth LL when we observe 
no features corresponding to scattering between nonequivalent nodes. This should gives rise 
to a smoother spatial dependence of the LDOS than for the energies corresponding to higher 
Landau levels. 

We focus on a localized impurity, but we expect our results to be quite similar for an 
extended Coulomb impurity. The most significant difference will be a reduction of the ratio 
between the intensity of the features away from the center and at the center. It would be 
interesting to see what happens if other types of disorder are considered which affect not only 
the electronic density, but also the hopping parameters in the neighborhood of the impurity. 
Unlike the potential disorder, the hopping disorder couples the nonequivalent valleys even in 
the zeroth LL, and one expects to observe the corresponding inter- nodal scattering features 
in the FTSTS spectra for all LLs. Therefore we propose to use the ratio between the intra- 
nodal and inter-nodal scattering features not only as a good indicator of the extension of 
the impurity potential, but also in the zeroth LL as an indicator of the form of the impurity 
potential (potential disorder versus hopping disorder). 

The relevant physical regime for the magnetic field is between 20 — 50T, and the dominant 
features we describe (for the second LL) are expected to arise for energies of the tip matching 
the energies of the LL's, i.e. 200 — 4007776^. It would be interesting to study what happens for 
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energies of the tip situated between two Landau levels in the presence of multiple impurities. 
Because of disorder, the LL's are expected to broaden, which translates spatially into a LL 
energy that fluctuates spatially. Thus, for some energies, one may have overlapping contribu- 
tions from various LL's, as well as more complicated (and not so neat) features arising in the 
FTSTS spectra. 

We have found that the FTSTS spectra can give information about the electronic wave- 
function in graphene at energies matching the LL energy (decay length, nodes, etc.), as well as 
about the nature of the disorder. We hope that a comparison between this theoretical study 
and experiments will shed light on the nature of the electronic states in epitaxial graphene 
under high magnetic field, where quantum-Hall features such as the LL's in the DOS are 
present, but no quantum Hall effect is observed. In particular it would be interesting to test 
whether the quasiparticles in epitaxial graphene are indeed described by typical quantum- 
Hall wavefunctions, and if so, whether it is possible to establish the nature of disorder in 
both epitaxial and in exfoliated graphene. Identifying the form of the wavefunction, as well 
as the nature of disorder will help understand why the quantum Hall effect is not observed in 
epitaxial graphene. 
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